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Abstract 

We complete here our work on isospin violation in the K ^ 3tt system. We first calculate K ^ 2n to the same 
order as we did K ^ 3tt in papers I and II of this series. This adds effects of order G27P^{iTiu — md) and G27P^e^ 
to earlier work. We calculate also the lowest order Bremsstrahlung contributions, K ^ttj, Sttj. With these and our 
earlier results we perform a full fit to all available CP conserving data in the K 27r, Stt system including isospin 
violation effects. We perform these fits under various input assumptions as well as test the factorization and the vector 
dominance model for the weak NLO low energy constants. 
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Abstract. We complete here our work on isospin violation in the K Sn system. We first calculate 
— !■ 27r to the same order as we did K ^ 3tt in papers I and II of this series. This adds effects of 
order G27 {rriu — md) and Girp^e^ to earlier work. We calculate also the lowest order Bremsstrahlung 
contributions, K — > 27r7, 3tt^. With these and our earlier results we perform a full fit to all available CP 
conserving data in the K 2t7, Stt system including isospin violation effects. We perform these fits under 
various input assumptions as well as test the factorization and the vector dominance model for the weak 
NLO low energy constants. 
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1 Introduction 

Low-energy QCD is non-perturbative, which calls for al- 
ternative methods of calculating processes including com- 
posite particles such as mesons and baryons. A method 
used to describe the interactions of the light pseudoscalar 
mesons (i^, tt, 77) is Chiral Perturbation Theory (ChPT). 
It was first presented by Weinberg, Gasser and Leutwyler 
[1-3] and it has been very successful. Pedagogical intro- 
ductions to ChPT can be found in [4]. The theory can be 
extended to also cover the weak interactions of the pseu- 
doscalars, first done in [5]. 

The first calculation of a kaon decaying into pions 
(K —^ 27r, Stt) was presented in [6], and reviews of other 
applications of ChPT to nonleptonic weak interactions can 
be found in [7]. 

The details from [6] were lost, but a recalculation in 
the isospin limit of K ^ 2tt to next-to-leading order was 
made in [8,9] and of if Stt in [9, 10]. In [9] a full fit 
to all experimental data existing at the time was made, 
and it was found that the decay rates and linear slopes 
agreed well. However, a small discrepancy was found in 
the quadratic slopes, and this can have several different 
origins. It could be an experimental problem or it could 
have a theoretical origin. In the latter case the corrections 
to the amplitude calculated in [9] are threefold: strong 
isospin breaking, electromagnetic (EM) isospin breaking 
or higher order corrections. 

In [11] the strong isospin and local electromagnetic 
corrections were investigated and it was found that the 

* Supported in part by the European Union TMR network, 
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inclusion of those led to changes of a few percent in the 
amplitudes. The local electromagnetic part was also calcu- 
lated in [10], in full agreement with our result after sort- 
ing out some misprints in [10], corrected in [12]. In [13] 
the radiative corrections were added, which means that 
the full effects of isospin breaking were studied. This led 
to changes in the amplitudes of order 5-10 percent. Note 
that the results in [13] disagree numerically with the re- 
suhs for K+ -> 7r07r07r+ of [14]. 

To answer the question whether isospin breaking re- 
moves the problem of fitting the quadratic slopes, a new 
full fit has to be done. That is the main result in this pa- 
per, and in this new fit we also include new experimental 
data [15, 16]. We also present recalculations of the ampli- 
tudes K 2tt, K — > 27r7 and K 37r7, all calculated 
to next-to-leading order and including first order isospin 
breaking, i.e. we include contributions proportional to p^, 
m^, e^, rriu — nid (leading order), and p'^,p^ m^, ^p^ e^, 
m^e^, p^{mu — nid) and ■m^{mu — md) (next-to-leading 
order). The corrections needed to be added to determine 
TTTT scattering lengths from K 3tt [17] are beyond the 
order calculated in this paper. 

The outline of this paper is as follows. The next section 
describes isospin breaking in more detail. In section the 
basis of ChPT, the Chiral Lagrangians, are discussed. Sec- 
tion^lspecifies the decays and describes the relevant kine- 
matics. The divergences appearing when including pho- 
tons are discussed in sectional In sectional the analytical 
results are discussed, section [7| contains the numerical re- 
sults and the last section contains the summary. 
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2 Isospin Breaking 

Isospin symmetry is the SU(2) symmetry under the ex- 
change of up- and down-quarks. This symmetry is only 
exact in the approximation that m„ = rrid and electro- 
magnetism is neglected, i.e. in the isospin limit. Calcula- 
tions are sometimes performed in the isospin limit since 
this is simpler and gives a good first estimate of the re- 
sult. However, to get a more accurate result one should 
include the effects from niu 7^ rrid and electromagnetism, 
i.e. isospin breaking. 

The two different sources of isospin breaking give rise 
to different effects. Strong isospin breaking, coming from 
TO„ 7^ md, include mixing between tt" and 77. This mix- 
ing leads to changes in the formulas for both the physical 
masses of tt" and 77 as well as the amplitude for any pro- 
cess involving either of the two. For a detailed discussion 
see [18]. 

The other source, electromagnetic isospin breaking, 
coming from the fact that the up- and the down-quarks are 
charged, implies interactions with photons. This means 
both the addition of new Lagrangians at each order, as 
well as the introduction of new diagrams including explicit 
photons. 



where the special unitary matrix U contains the Goldstone 
boson fields 



U — exp 



iV2 



M 



M 



V6 



V2 



K+ 



V 



K- 



(4) 



We use the formalism of the external field method [2] , and 
to include photons we set 



rrid 



and 



ieQAJJ - ie UQA^ 



where is the photon field and 



'2/3 



-1/3 



-1/3, 



(5) 



(6) 



(7) 



3 The ChPT Lagrangians 

The basis of our ChPT calculation is the various Chiral 
Lagrangians of relevant orders. We work to leading order 
in m„ — rud and but next-to-leading order in and 
m^. For simplicity we call in the remainder terms of order 
p 

p" TO', TO', p 

next-to-leading order 



TO^, e^, TO„ — m.d leading order and terms of order p'^, 
2 ^2 p2 g2^ ^2 g2^ p^irriu — md) and m?{mu — rud) 



The quadratic terms in |(2Jl are diagonalized by a rotation 



77 



7r3 cos e + r/g sm e 
— TTs sine -I- r/s cose , 



where the lowest order mixing angle e satisfies 
tan(2e) = Vi- 



md 



2 TO, 



md 



(8) 



(9) 



The weak, AS — 1, part of the Lagrangian has the 
form [19] 



3.1 Leading Order 

The lowest order Chiral Lagrangian is divided in three 
parts 

^2 = -^32 + Cw2 + ^E2, (1) 

where £52 refers to the strong AS = part, Cw2 the 
weak AS — ±1 part, and £e2 the strong-electromagnetic 
and weak-electromagnetic parts combined. For the strong 
part we have [2] 



Cs2 = ? {UnU'' 



^W2 = C Fq 



Gs{As2U^un + G'^{A32X+) 



+G27f^'" {A,,u^){Amu^) 
The tensor f^^'^^ has as nonzero components 



^21,13 _ ^13. 21 



^23,33 _ ^33,23 



_|_ j.22,23 _ i23.22 

3 ' 

}_ ^23,11 _ ^11'23 
6 ' 



h.C. . 



1 

^6 

1 

3 ' 



(10) 



X4 



and the matrix An is defined as 



Here (A) stands for the flavour trace of the matrix A, and 
Fq is the pion decay constant in the chiral limit. We define 
the matrices u^, u and x± 



(3) 



Aij = uXijv) , {■^ij)ab ~ ■ 



(11) 



The coefficient C is defined such that in the chiral and 
large Nc hmits Gg = G27 = 1, 



C 



3 Gp 

"5 71 



Vud V,* 



-1.06- 10"'' GeV" 



(12) 
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Finally, the remaining electromagnetic part, relevant for 
this calculation, looks like (see e.g. [20]) 

£b2 = e^F^Z{QLQR) + e^F^iTQu) , (13) 

where the weak-electromagnetic term is multiplied by a 
constant Ge (ffewkGg in [20]), 

T = Ge F^A32 + h.c. (14) 

and 

Ql = uQuU Qr = v)Qu. (15) 
3.2 Next-to-leading Order 

Chiral Perturbation Theory is a non-renormalizable the- 
ory. This means that new terms have to be added at each 
order to compensate for the divergences coming from loop- 
diagrams. Thus the Lagrangians increase in size for every 
new order and the number of free parameters rises as well. 
At next-to-leading order the Lagrangian is split in four 
parts which, in obvious notation, are 

i^S4,+ i^W4:-\- ^S2E2+ i^W2E2{Gs) ■ (16) 

The notation (Gg) indicates that here only the dominant 
Gg-part is included in the Lagrangian and therefore in the 
calculation. 

The Lagrangians of next-to-leading order are quite large 
and we will not write them explicitly here since they can 
be found in many places [2,5,20-24]. For a list of all the 
pieces relevant for this specific calculation see [11,13]. Note 
however that one contributing term was forgotten when 
writing Cs2E2 in [11], namely 

-V^,QrQr + QrV^,Qr)u''), (17) 

where 

V^Ql = V^Qi + ^[u^, Ql] = uD^Qlv), 
V^Qii = V^Qij-|[w^,Qij]=«tD^Qiju. (18) 

It contributes to the calculation of the decay constants, 
F^+ and Fji+. It only contributes to the amplitudes of 
K ^ 27T and K ^ 3tt via the rewriting of the lowest 
order in terms of F^+ and Fx+ rather than Fq. 

3.2.1 Ultraviolet Divergences 

The processes K ^ 27r and K Sir receives higher-order 
contributions from diagrams that contain loops. The study 
of these diagrams is complicated by the fact that they 
need to be precisely defined. The loop-diagrams involve an 
integration over the loop- momentum Q, and the integrals 
are divergent in the ultraviolet region, i.e. when Q — > oo. 
These ultraviolet divergences are canceled by replacing the 
coefficients, Xi, in the next-to- leading order Lagrangians 
by the renormalized coefficients, XI, and a subtraction 
part. See [9,11] and references therein. 



4 Kinematics 

4.1 K an6 K ^ 2-K'y 

In the limit of CP-conservation, there are three different 

decays of the type K ^ 2tt {K~ decays arc not treated 
separately since they are counterparts to the decays): 

X+(A:)^^+(pi)^"(p2), [A+o], (19) 

where we have indicated the four-momentum defined for 
each particle and the symbol used for the amplitude. With 
an external photon it changes to: 

Ks{k)^7r°{p^)7r°{p2)l{q), [AI^,] , 
Ks{k)^n+{p,)7T-{p2)j{q), [Al_^], 
K+{k)^n+{p,)n"(p2)j{q), [A+o^] . (20) 

The kinematics for K 27r7 is treated using 

ro = -k-q, ri=pi-q, r2=P2-q, (21) 

where 

ro + ri+r2=0. (22) 

4.2 K ^Stt and K ^ Swj 

For the corresponding process K 37r, there are five 
different decays: 



KLik) - 


-.7r0(pi)7r0(p2)7r0(p3), 


[^000 ] ' 




KL{k) - 


7r+(pi)7r"(p2)7r°(p3) , 






Ks{k) - 


7r+(j?i)7r"(j?2)7r°(p3) , 






K+{k) - 


-7r0(pi)7r0(p2)7r+(p3), 


[^00+] , 




K+{k) - 


7r+(]3i)7r+(p2)7r"(p3) , 


[A++-], 


(23) 



and here the variables are 

Sx = {k-pxf , S2=(k-P2f , S3 = {k-P3f, 

(24) 

where 

si + S2 + S3 = k^+pl+pl+pl. (25) 

The amplitudes are expanded in terms of the Dalitz plot 
variables x and y defined as 

S3 - So S2 - Si 1 
y = 2 ' ^ = 2 ' So = - (Si -h S2 + S3) . 

m + m + 3 

(26) 



With an external photon the decays are: 




KUk) - 


■.7rO(pi)7rO(p2)7rO(p3)7(5), 


[^0007] ' 


KL{k) - 


7r+(pi) 7r"(p2) 7r°(p3) 7 {q) , 


[^+-07] > 


Ks{k) - 


7r+(pi)7r"(p2)7r°(p3)7(g) , 


[^+-07] 5 


K+{k) - 


-7r"(pi)7r"(p2)^+(p3)7(9), 


[A0O + 7] , 


K+{k) - 


7r+(pi)7r+(p2)7r"(p3)7(9) , 


[^++-7] .(27) 
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where the kinematics is treated using 

si^ = {k- pif , S27 = (fc - , = (fc - psf , (28) 

to = ~k-q, ti=pi-q, t2=P2-q, h=Pz-q, (29) 
where 

+ + t2 + t3 = (30) 

and 

S17 + S27 + S37 = fc^ + Pi + P2 + P3 ~ 2 io ■ (31) 



6 



5 Infrared Divergences 

Beside the uhraviolet divergences which are removed by 
renormahzation of the higher order coefficients, diagrams 
including photons in the loops contain infrared (IR) diver- 
gences. These infinities come from the Q ^ end of the 
loop- momentum integrals. They are handled by including 
also the Bremsstrahlung process, where a real photon is 
radiated off one of the charged mesons. It is only the sum 
of the virtual loop corrections and the real Bremsstrahlung 
which is physically significant and thus needs to be well 
defined. 

We regulate the IR divergences in both the virtual 
photon loops and the real emission with a photon mass 
ni'y and keep only the singular terms plus those that do 



not vanish in the limit 



0. We include the real 



Bremsstrahlung for photon energies up to a cut-off lu and 
treat it in the soft photon approximation. 

The exact form of the amplitude squared for the Brems- 
strahlung depends on which specific amplitude that is be- 
ing calculated. For a detailed presentation of the calcula- 
tion and resulting expressions for if — > 37r see [13]. The 
corresponding amplitudes for K 2n are 

0, (32) 



-IBS 



-\LO 



47r2 



log - IiR{m^,m^, tjik) 
mz 



(33) 



l^+o|ls 



\A+o\' 



LO 



47r2 



log- 



r / 2 2 2\ 



where 

IiR{mj,ml,ml) = 



(34) 



Xs m\ — TO? — Too , , up 

71 rloga^s log— 5-. 

47r"^ TOiTO2(l — Xs) mi^ 

(35) 

When using these expressions, the divergences from the 
photon loops cancel exactly. 

A similar problem shows up in the definition of the 
decay constants, since we normalize the lowest order with 
J^7r+ and Fk+ . See [13] for details. 



6 Analytical Results 

6.1 K ->27r 

The most complete work on isospin violation in K —t 2tt 
is in [25], earlier work can be found in [26]. 



Fig. 1. The diagrams for K —* 2-n. An open square is a ver- 
tex from Cw4, or Lw2E2, a filled square a vertex from Cw2 
or Ce2 {^S — 1) and a filled circle a vertex from £52 or 
jC-e2 {^S = 0). a straight line is a pseudoscalar meson and 
a wiggly line a photon. 



6.1.1 Lowest Order 

There is only one diagram contributing to the decay K 
2tt at lowest order, see top left in Fig.^ and the resulting 
amplitudes are also quite simple. To first order in isospin 
they can be written 



A^O = ZF(G8-G27) 

/ . sine 



(36) 



V3 



(TO^-TO;^)-2TO^o+2TO^a , (37) 



iFGs (-2to2+ +2to|-o) 



+lFG27 



-3 + 3 ™ 



K+ 



2iF^e^ Ge 



(38) 



sme , „ sme 



A+o^iFGs { -2 ^ ml + 2 ^ m\ - -|- ml„ 



I ■ j?n \ Q 2 , Q sine 2 7/02 

+iFG2T I - 3 TO^ -t- 3 TO^ - 7/3 

+2/3mlo+5/3ml^^ -iF^e'^GE- (39) 
See section 17. 1.1 1 for a discussion of the masses used. 



6.1.2 Next-to-leading Order 

There are seven more diagrams contributing to next-to- 
leading order, see Fig. ^ The resulting amplitudes are 
long, and we decided to not include them here but in- 
stead make them available for download [27]. Note that 
we have also included contributions proportional to G27, 
not included in [25]. These are included for consistency 
between the K 27r(7) and K 37r(7) calculations. 
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K 



Fig. 2. The diagrams for K — » 2tt^. A square is a weak vertex, 
a circle a strong vertex, a straight line a pseudoscalar meson 
and a wiggly line a photon. 



6.2 K -> 27r7 

The amplitudes for the processes K — > 2tt'^ have been cal- 
culated before. Here we only need the lowest order contri- 
bution to be consistent with the K ^ 2tt calculation. This 
we recalculated and the starting point is the two diagrams 
that contribute to the process, shown in Fig. [3 

Since photons at this order only couple to charged par- 
ticles, there are only two different amplitudes and the re- 
sults are 



— Got 

3 



and 









-a 






6 


(to|- - 




1 


+Pi-e 






ro 



{m\ - ml) 



Pi.e) 



it 









(40) 



—k.e 



P2-S 



1 

ro 
1 



1 
1 



(.41) 



These amplitudes can be decomposed into an electric and 
a magnetic part: 



A{K 27r7) = ee''{k) [E^ 



T\/rvp<y\ 



(42) 



but at lowest order the magnetic amplitude M'^'"^ vanishes 
since there is no Sf^iipa tensor in the corresponding lowest 
order Lagrangian. 

The electric amplitude, on the other hand, is com- 
pletely determined by the corresponding non-radiative am- 
plitude via Low's theorem [28]. 



6.3 K ^ 37r7 

The decay K 37r7 is discussed in detail in [29] . We only 
need the lowest order amplitude for consistency with the 
calculation oi K ^ 3tt. We have calculated the four dif- 
ferent amplitudes using Chiral Perturbation Theory, and 
checked that they agree with Low's theorem. The calcu- 
lation is based on seven diagrams, see Fig. |2| The four 
amplitudes are 



A 



4—07 



Gg — G27 




Fig. 3. The diagrams for K Sny. A square is a weak vertex, 
a circle a strong vertex, a straight line a pseudoscalar meson 
and a wiggly line a photon. 



Gs G27 -3mj^ + 8ml 



3 6 

(~3(s37-2<i-2<2)+m^ + 3m^) x 
pi-e P2-£^ 



h t2 

Gs G27 -3 m% + 8 ml 
~3 



it2+tl) 



6 



Pi-e p2- e 



tl 



t2 



(43) 



-07 



G 



= I e — — 
6 



+2 to 
+2 to 



27 



m^ — 




S27) 


'Pl ■ £ 

. tl 


Pi ■ £ 


k ■ e 


tl 




P2 ■ £ 


k ■ £ 


t2 





{2 ml - 3to|-) X 



(44) 



A 



00+7 



^8.2, 2\ ^27 

2 ml — n 



K 



-13/6mt-M/2m],)-(% + 



(5/3 m\ 
G27 



mt — m 



K 



(7/6m2 - 17/6m^)) x 
(2 (S3^ - 2ti - 2t2) - m\ -iml) 
— 5/6(2 m^- 3 ml-) x 
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(-(S3^ - 2 ti - 2 ta) + mj^ + 3 ml) 
k ■ e ^ P3 ■ £ 



(45) 



Gg (ml + m\) + G27 {m\ + 13/3 m^) 



-(G8-13/3G27)(s37- 2^1-2^2) 



Ge 


-0.4 




0.38 


10" 


-3 


sin e 


1.19 ■ 10"^ 


i2 


1.59 


10- 


-3 


Z 


0.805 


LI 


-2.91 


10- 


-3 




0.0924 GeV 


LI 









Fk 


0.113 GeV 


LI 


1.46 


10- 


-3 




-10.4 • 10"^ 


LI 











5.95 • 10"^ 




-0.49 


10- 


-3 


Dl3 





LI 


1.0 


10- 


-3 


Dl5 





LI 


7.0 


10- 


-3 



Table 1. The various input values used, LECs given at 
0.77 GeV. 



k ■ £ ^ Pi ■ e ^ P2 ■ e 

to ti t2 



P3 ■ £ 



+(2G8-26/3 G27) [h-h) 



Pi-e p2-e 



t2 



(46) 



Once again, these amplitudes can be decomposed into an 
electric and a magnetic part: 



A{K ^ 3777) = e s^{k) {E^ + e^^p^^rP") , 



(47) 



and the magnetic amplitude M^P'^ vanishes for the reasons 
given above. 

The electric amplitude at this order is again completely 
determined by the corresponding non-radiative amplitude 
A{s, v) via Low's theorem [28,29]: 



EP = A{s,v) UP 

ds 
+ 0(k) 



A^2 + '-^iA^.. 



dv 



(48) 



with (the meson charges in units of e are denoted qi , with 



s 

V 

SP 



{pi +P2) 
k- {pi~ P2) 
4 



i=l 
ji 



AP, = {q,t, - q,ti)DP^ 



u 



(49) 



Since there are no terms of 0{k) at lowest order in the 
chiral expansion, the leading-order electric amplitude is 
completely determined by the explicit terms in H48|l . 



7 Numerical Results 

7.1 Experimental Data and Input 

For the numerical studies we use the input given in Ta- 
ble [3 



7.1.1 Strong and Electromagnetic Input 

There are different ways to treat the masses, especially in 
the isospin limit case. In [9] the masses used in the phase 
space were obtained from the physical masses occurring in 
the decays. However, in the amplitudes the physical mass 
of the kaon involved in the process was used and the pion 



mass was given by = 



=1,3 ' 



with i — 1, 2, 3 



being the three pions participating in the reaction. This 
allowed for the correct kinematical relation si + S2 + S3 = 
+ 3m^ to be satisfied while having the isospin limit in 
the amplitude but the physical masses in the phase space. 
The results in [9] were obtained with the physical mass 
for the eta. Results with the Gell-Mann-Okubo (GMO) 
relation for the eta mass in the loops gave small changes 
within the general errors given in [9]. 

In the decays here, we work to first order in isospin 
breaking. We have rewritten explicit factors of m„ — 
in terms of sin e according to 



1 



ma 



md) sme . 



(50) 



In general we use the physical masses of pions and kaons 
in the loops but as soon as a factor of sin e or is present 
we use a common kaon and a common pion mass. This 
simplifies the analytical formulae enormously. The kaon 
mass chosen is the mass from the kaon in the decay and 
the pion mass used is 3m^ = m^; with i = 1,2,3 the 
three pions in the final state, i.e. the mass we used in 
the isospin limit case. For the eta mass we use in general 
the physical mass in the loop integrals. We have used the 
GMO mass relation with isospin violation included, 

2 1 

= g ("1a'+ + ~ ) + 3 "^7r« > (51) 

to simplify the amplitudes, except in the loops as stated 
above. The possible lowest order contributions from the 
eta mass have been removed from the amplitudes using the 
corresponding next-to-leading order relation as described 
in [11]. 

The strong LECs, L\ to Lg, as well as sine come from 
the one-loop fit in [18], from [30] and the Ge estimate 
is from [31]. 

The constant Z from Ce2 we estimate via 



Z = 



1 



2F2, 



("4+ - "^^0), 



(52) 
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which corresponds to the value in Table ^ The higher 
order coefficients oi Cea, . . . are rather unknown. 
Some rough estimates exist but we put them to zero here, 
at the relevant scale. 

The IR divergences are canceled by adding the soft- 
photon Bremsstrahlung. We have used a 1 MeV cut-off 
in energy for this and used the same cut-off in the defini- 
tion of and Fx+- We also use rrij = 1 MeV, which 
effectively removes the infrared part. 

The subtraction scale /i is chosen to be 0.77 GeV unless 
stated otherwise. 



7.1.2 Input Relevant for the Photon Reducible Diagrams 

The two constants D13 and D15 are set to zero since no 
knowledge exist of their values. One can determine the 
constants A^i4 and A^i5 from K Trl~^l~ decays. For a de- 
tailed analysis, see [13]. The resulting values are given in 
Tabled Note however, that as described in [13], D13, D15, 
Ni4 and Ni^ only contributes via the photon reducible 
diagrams. These diagrams are negligible numerically, un- 
less the constants are orders of magnitude larger than ex- 
pected. 



7.2 Bremsstrahlung and Dependence on and oj 

The isospin breaking amplitudes for K 2tt and K —^ Sir 
both depend on , introduced to regularize the infrared 
divergences coming from loops containing photons. This 



Lo (GeV) 


IR photon 


Extra photon 


Sum (GeV) 




(GeV) 


(GeV) 




























0.01 


4.34 


• 10"^' 


1.59 


10" 


iV 


5.92 


10" 


17 


0.001 
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3.69 


10" 


17 


5.86 


10" 


17 


0.0005 


1.52 




4.34 


10" 


17 


5.85 


10" 


17 


0.0001 







5.85 


10" 


17 


5.85 


10" 


17 






K+ 












0.01 


4 'Hn 


■ 10-^" 


1.54 


1 A~ 


•2U 


5.84 


10" 


20 


0.001 


2.15 


• 10-2" 


3.61 


iU 


20 


5.77 


10" 


20 


0.0005 


1 


■ 10-2" 


4.26 


1 n~ 

±u 


20 


5.76 


10" 


20 


0.0001 







5.77 




20 


5.77 


10" 


20 






Kl - 




TV 










0.01 




■ 10-'^' 


3.39 


1 A — 
10 


22 


2.73 


10" 


21 


0.001 


1 1 Q 


■ 10-21 


1.39 




21 


2.58 


10" 


21 


0.0005 


R 'H4 


• 10-22 


1.74 


10" 


21 


2.58 


10" 


21 


0.0001 







2.56 


10" 


21 


2.56 


10" 


21 


















0.01 


0.00 


■ 10-2^ 


8.93 


10" 


2!) 


6.77 


10" 


24 


0.001 


2.94 


• 10-24 


3.43 


10" 


24 


6.37 


10" 


21 


0.0005 


2.05 


• 10-24 


4.28 


10" 


24 


6.34 


10" 


21 


0.0001 







6.29 


10" 


24 


6.29 


10" 


24 






K+ - 


— ^ TT TT 


n+ 










0.01 


3.77 


■ 10-^^ 


4.03 


10" 


23 


4.18 


10" 


22 


0.001 


1.89 


• 10-22 


1.98 


10" 


22 


3.86 


10" 


22 


0.0005 


1.32 


• 10-22 


2.53 


10" 


22 


3.85 


10" 


22 


0.0001 







3.83 


10" 


22 


3.83 


10" 


22 






K+ - 














0.01 


5.05 


■ 10-21 


5.86 


10" 


22 


5.63 


10" 


21 


0.001 


2.52 


• 10-21 


2.73 


10" 


21 


5.25 


10" 


21 


0.0005 


1.76 


• 10-21 


3.46 


10" 


21 


5.23 


10" 


21 


0.0001 







5.19 


10" 


21 


5.19 


10" 


21 



-dependence is canceled by adding the bremsstrahlung Table 2. K ^ 2, 37r decay rates calculated for different values 



amplitudes, where a real soft photon is radiated off one 
of the mesons. This cancellation was checked for K ^ Sir 
in [13], and we have now also checked it for K — + 2tt. 

However, after the addition of bremsstrahlung the de- 
cay rates depend instead on w, the cut-off energy of the 
radiated real photon. This is a parameter that should be 
set to a value depending on the experiment that one com- 
pares to. 

Another possibility, which we use in this paper, is to 
add the full amplitudes with an extra radiated photon, 
K 27r7 and K Sttj. When doing that the decay 
rates should be independent of uj. We have checked this 
numerically and the results are presented in Table |5J For 
this comparison we have chosen m-y = 0.1 MeV and varied 
omega over a large range. One can see that up to photon 
energies of 1 MeV, the sum is constant within the expected 
uncertainties. The different sum when including energies 
up to 10 MeV is an indication that the soft photon approx- 
imation, used in calculating the infrared contribution, is 
breaking down. 

The way we treated the Bremsstrahlung contribution 
in the fits is as follows. We assume that the measured de- 
cay widths are including all photons. To compare to our 
amplitudes (calculated without hard photons), we there- 
fore subtract numerically the calculated hard photon con- 
tributions from the experimental numbers. 



on LO. Here we use m-y 



0.1 MeV and the same value is used 



as a cut-off in the decay constants, see the appendix in [13] 



7.3 Fit to K 



2n 



The process X — > 27r in the presence of isospin breaking 
has been discussed in detail in [25]. We have reproduced 
that calculation but added in addition also all the isospin 
breaking contributions from the 27 amplitudes, except for 
the parts from the weak-electromagnetic 27 Lagrangian. 

The isospin breaking corrections to the decay rates are 
rather small, but they have an impact on the phase shift 
between the 1 = 2 and / = amplitudes, S2 — Sq, as 
described in detail in [25]. Our results are compatible with 
the ones presented there. The phase shift we use here is 
defined via 



s V2 2 



A 



V3 
V3 



V3 



-1-0 
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K ^2Tr 


Order 


/i [GeV] 


Gs 


G27 


S2 — So 


LO 




10.4 


0.55 


-60.3° 


NLO 


0.6 


6.43 


0.44 


-58.9° 


NLO 


0.77 


5.39 


0.36 


-58.0° 


NLO 


1.0 


4.60 


0.30 


-57.4° 



Table 3. Quantities fitted from K 2n only. Tlie scale /j. is 
the scale at which the unknown low energy constants are put 
to zero. 



A2 
Ao 



^0 



(53) 



In the fits we have left 82 — 5a as an additional free param- 
eter, as it is known that this phase is badly reproduced at 
one-loop order in ChPT. Note that because of the isospin 
breaking the A2 amplitude is different from A2. 

We have performed a lowest order and a NLO fit to 
only the K ^ 2ti amplitudes. In the NLO fit we set all 
the extra parameters, Ki = = Zl = at the scale 
jj, indicated. The Bremsstrahlung contribution has been 
subtracted as discussed above. 

As can be seen in Table O there is a sizable vari- 
ation depending on the input scale used. There is very 
little change in the absolute values of Gg and G27 com- 
pared to the isospin limit fit of [9] , where only the fit with 
1^ = 0.77 GeV was done. The angle is similar to the fit 
there, but this is a combination of two different effects. It 
was lowered because the new KLOE data have now been 
included in the PDG averaging, but the isospin break- 
ing effects induced a positive correction as was also found 
in [25]. 

The values of Gs and G27 are determined by fitting 
CF^^Gs and CF^G27 and then setting Fq = F^ numeri- 
cally to provide the numbers in the tables. 



7.4 Fit to K ^2t: and K 



The quantities we fit are the measured decay rates and the 
various parameters of the Dalitz plot distributions defined 
via 

^(si, 52,53) 



A{so,so, So) 



1 + gy + hy + kx 



(54) 



For the decay Kl Stt^, k = h/3 and g = 0. The decay 
Ks 7r"'"7r~7r'^ is included via 



A 



s 



= isx - ixy . 



(55) 



The decay rates are included in the fit as follows. We 
subtract from the decay rates the Bremsstrahlung contri- 
butions as described above as a function of Gs and G27. 
We then convert the decay rate using the central values of 
the measured Dalitz plot distribution into a value for the 
amplitude squared at the center of the Dalitz plot. These 
squared amplitudes together with the parameters g, h, k 
and 7s are used as the 18 experimental parameters to be 
fitted. 



This means that the effect of Bremsstrahlung is in- 
cluded fully in the decay rates, but only via the soft pho- 
ton approximation with a 1 MeV cut-off for the Dalitz 
plot distributions. We have not included the preliminary 
data from KTeV, NA48 and KLOE. 

The number of free input parameters on the theory side 
is very large. Since it turns out that the isospin breaking 
effects are very small, we put those extra NLO parameters 
equal to zero at the scale /i indicated in the tables. 

Let us repeat here the definitions of the various ex- 
tra NLO parameters. The i[ are taken from the stan- 
dard fit done at one loop to be compatible with the order 
of this calculation. The are the extra parameters at 
NLO in the p^e^ sector. Those are always put to zero at 
the scale indicated. In the isospin limit 11 combinations 
of the weak NLO low-energy coefficients show up, as dis- 
cussed in [9]. These are Ki, i — 1, . . . , 11. In the presence 
of isospin breaking many more combinations of these, as 
well as from the weak octet order e^p^ Lagrangian, emerge 
and they were classified in [11]. The 27-part of the weak 
Lagrangian of order e'^p^ has not been worked out and will 
lead to some more free parameters. We have not used any 
estimates of these extra parameters but set all of them 
to zero at the scale indicated, except for Ki, i — 1, . . . , 7. 
The reason for this choice is that they are the leading con- 
tributions. i^i,2,3 are octet enhanced and come multiplied 
with factors of order and -ftr4,5,6,7 are 27-plets but also 
come multiplied with factors of order . The neglected 
ones are thus suppressed by either isospin breaking, fac- 
tors of m^/m|. or by the AI — 1/2 rule, i.e. an extra 
factor of G27/G8- 



7.4.1 General Fits 

Here we perform the fits with similar assumptions as used 
in the isospin limit fit, as well as a few additional ones. 
First Gs and G27 are extremely correlated with the values 
of Ki and K4 respectively. They are very difficult to obtain 
separately without additional assumptions. The main fit 
is therefore the one with 



Ki^Ki^Ks^Kg^O, 



(56) 



at a scale ^ — 0.77 GeV. The results are given in TableEl 
This table is very similar to Table 6 in [9] . The large values 
of Kq and the resulting large value of K7 have the same 
origin as in that reference. In order to fit 75 well, Kq is put 
large because it gets multiplied there with a small factor 
and is the only one contributing. This in turn leads large 
values for K7 to compensate in other places. 
The fit with 

K6 = (57) 

in addition has only a slightly larger and a smaller 
K7. The is larger than in [9] because the experimental 
errors on several quantities have decreased since then. The 
overall fit is slightly better than the one of [9] because the 
newer measurements of the Dalitz distribution in K^ —>■ 
TT^Tr^Tr"*" agree better with the chiral expressions. 



8 Johan Bijnens, Fredrik Borg: Isospin Breaking in K ^ 3n Decays III: Bremsstrahlung and Fit to Experiment 



Decay 


Width [GeV] 


ChPT [GeV] 


Fact. [GeV] 


A5 — TT TT 

Ks tt^tt" 


(1.1231 ± 0.0078) ■ 10"'' 
(2.2828 ± 0.0104) ■ 10"^^ 
(5.0691 ± 0.0108) -10"^^ 


1.123 • 10"'' 
2.282 • 10""^ 
5.069 ■ 10""^ 


1.127 ■ 10""' 
2.283 ■ lO""'' 
5.069 ■ lO-"'^ 




(2.6748 ± 0.0358) ■ lO"'* 
(1.5998 ± 0.0271) ■ 10"^** 
(9.195 ±0.0213) • 10""^ 
(2.9737 ± 0.0174) ■ 10"^** 


2.618 • 10"'** 
1.658 ■ 10""** 
8.934 • 10""^ 
2.971 • 10""** 


2.698 ■ lO"'** 
1.711 ■ 10""* 
8.816 ■ 10""^ 
2.933 ■ 10""* 



Table 4. The various decay widths from the PDG tables [32], and our results from the main fit and the best factorization fit. 



Decay 


Quantity 


Experiment 


ChPT 


Fact. 


Kl 


— > TT TT TT 


h 


-0.0050 ± 0.0014 


-0.0062 


-0.0025 


Kl 




g 


0.678 ± 0.008 


0.678 


0.654 






h 


0.076 ± 0.006 


0.088 


0.083 






k 


0.0099 ± 0.0015 


0.0057 


0.0068 


Ks 




-ys 


(3.3 ±0.5) ■ 10"** 


3.0 ■ 10"" 


2.9 • 10"** 


K^ 


— > TV TT TV 


g 


0.638 ± 0.020 


0.636 


0.648 






h 


0.051 ±0.013 


0.077 


0.080 






k 


0.004 ± 0.007 


0.0047 


0.0069 


K+ 


TT^TV~^TV~ 


g 


-0.2154 ± 0.0035 


-0.215 


-0.226 






h 


0.012 ± 0.008 


0.012 


0.019 






k 


-0.0101 ± 0.0034 


-0.0034 


-0.0033 


K- 


7r"7r"7r^ 


g 
h 
k 


-0.217 ±0.007 
0.010 ± 0.006 
-0.0084 ±0.0019 







Table 5. Experimental values and the main fit and best factorization fit of the Dalitz plot distribution parameters. The data 
are from the PDG tables [32] except 7s from [33]. 



Constraint 


Eq. m 


Eq. m 


Eq. m 


Eq. m^il 




0.77 GeV 


1.0 GeV 


0.6 GeV 


0.77 GeV 


Gg 


5.39(1) 


4.60(1) 


6.43(1) 


5.39(1) 


G27 


0.359(2) 


0.301(1) 


0.438(2) 


0.359(2) 


S2 — So 


-57.(1.5)° 


-57.3(1.4)° 


-58.9(1.4)° 


-57.9(1.4)° 


W^'Ki/Gs 


= 


= 


= 


= 


10^7^2 /Gg 


48.5(2.4) 


56.5(2.4) 


41.2(1.9) 


46.6(1.6) 


lO^T^a/Gg 


2.6(1.2) 


-1.7(1.1) 


6.7(1.0) 


3.5(0.8) 


10^7^4/G27 


= 


= 


EE 


EE 


m^Ks/G27 


-41.2(16.9) 


-52.0(17.7) 


-31.1(12.0) 


-27.0(8.3) 


10'K(i/G27 


-102(105) 


-114(105) 


-93(76) 


= 


I(fk7/G27 


78.6(33) 


78.0(33.5) 


79.6(22.7) 


50.0(13.0) 


lO^A-g/Gg 


= 


= 


= 


= 


IQ^kg/Gs 


EE 


= 


= 


EE 


X7D0F 


29.3/10 


27.2/10 


33.0/10 


30.5/11 



Table 6. The results for Gg and G27 and the Ki for the various constraints described in the text. In brackets are the MINUIT 
errors. The Ki are quoted at the scale ^ mentioned. 



We get fits of roughly similar quality for all values of 
jjL where the other parameters have been put to zero. The 
fits tend to be slightly better for the larger values of /i. 
The fitted values for the Ki are /^-dependent, albeit not 
extremely strongly. 

The Ki themselves have a /i-dependence which is given 
by the cancellation of divergences, and this can be calcu- 
lated from the known subtractions. We have shown the 



variation with /i from /i = 0.77 GeV to fi = 0.6 GeV and 
/I = 1.0 GeV for K,, i ^ 1, . . . ,11 in Table[3 

In order to compare with the factorization model of the 
weak low energy constants, we also perform a fit where all 
next-to-leading order LECs proportional to G27 are set to 
zero, but we keep in addition the sub-leading octet ones. 
This fit is shown for ^ ^ 0.77 GeV in Table [7| The fit 
is somewhat worse than those of Tab. but not much. A 
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Constraint 


Eq. m 


ji variation 


ji variation 


M 


0.77 GeV 


1.0 GeV 


0.6 GeV 


Gs 


4.84(1) 


- 


- 


G27 


0.430(1) 


- 


- 


S2 — So 


-57.9(0.2)° 


- 


- 


lO^ATi/Gg 


2.0(1) 


-5.88 


5.61 


10=^^-2 /Gs 


63.0(1.5) 


-2.69 


2.57 




-6.0(7) 


0.159 


-0.152 


W^K4/G27 


= 


-9.93 


9.48 




= 








lVpk6/G27 


= 


27.0 


-25.8 


10=^i^7/G27 


= 


-21.5 


20.5 


lO^^A-g/Gg 


20.4(1) 


-0.546 


0.521 


10=^7^9/G8 


9.1(1) 


-2.92 


2.79 


lO^ifio/Gg 


= 


11.6 


-11.1 


10='A-ii/G8 


= 


-1.66 


1.58 


x7DOF 


33.3/10 







Table 7. The results for Gs and G27 and the Ki for the octet 
constraint described in the text. In brackets are the MINUIT 
errors. The Ki are quoted at the scale ^ mentioned. The last 
two columns are the values of the Ki at the scale ^ mentioned 
when they are zero at = 0.77 GeV and run with Gs = 5.39 
and G27 = 0.359. 



very similar fit is obtained for n — 1 GeV with a of 
29.9. At /i = 0.6 GeV the best solution found had a of 
57.8. This fit corresponds to 

ki = kr, =kfi=I<7 = Q, (58) 

and this type of fit is referred to below as an octet fit. 



7.4.2 Fits to Factorization and other IVlodels 

Various models for the NLO weak constants exist. We will 
discuss here some of the ones which are presented in [23] . 
These have been discussed in that reference only for the 
pure octet case. So the quality of the models should be 
compared with the octet fit above. 

A first choice is the resonance exchange domination 
of the weak constants. The problem here is that the weak 
decays of the resonances involve themselves many new un- 
measured parameters and thus leads to fairly few general 
conclusions. If we assume that the vector octet exchange 
dominates, we get a relation between the octet NLO con- 
stants 

Nl + Nl^ +2Nl^{), (59) 

which is a combination we can in fact determine. It trans- 
lates for our parameters into 

kz = -^I<2. (60) 

It can be easily seen from Tables |B1 and [7] that this is very 
far from being satisfied by our fits. 

A very often used model is the factorization model. It 
corresponds to taking the underlying four quark operator 





0.77 GeV 


0.9 GeV 


0.842 GeV 


Gg 


4.18(1) 


4.42 


4.22(1) 


G27 


0.360(2) 


0.326(10) 


0.339(10) 


kp 


2.61(1) 


4.94(2) 


3.60(5) 


X7D0F 


109/14 


182/14 


60.4/13 



Table 8. The results for the fit with the factorization assump- 
tion for various values of n including the optimal one. 

and bosonizing separately the two quark currents present 
there. Looking at the dominant octet operators only for 
the cases that we need here, this leads to the relations [23] 

N{ = 2fc/(32/3i5; + 4L;; + 2/3L5), 

Nl^ = 2fc/(16/3i5; + ALl + 10/3 i§) , 

Nl = 2kf{SLl-2Ll), 

Nl = 2fc/(-16/3LJ-8/3i?;-4/3Lg), 

N;^2kf{~Ll), 

N; = 2kf{2/iLl^, 

m = 2kf{Ll), 

Nl = 2kf{^Ll + 2Ll), 

- N^o = Nil = = N^s = . (61) 

The parameter kf allows for some overall adjustment. The 
special value fc/ — 1/2 is referred to as the Weak Defor- 
mation Model (WDM) [23,34]. We have performed a fit 
leaving both kj and fi free. Note that the scale fj. is also 
the scale where we have put all the other NLO parameters 
equal to zero. The input values of the L[ have been scaled 
accordingly. 

The fits done with kf,Gs and G27 as free parameters 
have significantly larger than those reported above. 
Some representative values are shown in Table |H1 

The fit with /i free gave a minimum a.t n — 0.842 GeV. 
The fits with outside the range of Table|Hlhad very large 
values of x^- 

In order to show the quality of the fits we have given in 
Tables 0] and El also the values obtained for the quantities 
from the main fit and best factorization fit, labeled ChPT 
and Fact, respectively. Notice that the extrapolation to 
the full phase space here has been done from the ampli- 
tude squared in the center of the Dalitz plot using the 
experimental values for the distribution over the Dalitz 
plot. 

8 Summary 

We have recalculated in this paper the Bremsstrahlung 
amplitudes for K —^ 2TTJ and K — > 3777. In addition we 
have calculated also the isospin violating effects to K ^ 
27T including those with the 27-operators both for effects 
due to niu — ma and electromagnetism. This we did to be 
consistent with the calculations oiK ^ iir done in [11,13]. 

We have checked explicitly that the infrared diver- 
gences of the photon loops regulated by cancel between 
the virtual photon loops and the soft Brems-strahlung. 
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We checked in addition that the photon energy cut-off de- 
pendence cancels between the soft-photon Bremsstrahlung 
part and the part where hard photons are treated explic- 
itly. We have not included the explicit expressions for the 
K ^ 2tt amplitudes because they are rather long. They 
can be obtained from the authors or [27]. These ampli- 
tudes have passed all the standard tests, like cancellation 
of divergences from both NLO ChPT as well as the in- 
frared singularities. 

With these calculations and those published earlier 
in [11, 13], we have updated the fit to the CP conserv- 
ing observables in the K — > 7r7r(7r)(7) system done in 
the isospin limit in [9]. As expected from the fairly small 
isospin violating effects found in [11,13] and from the anal- 
ysis of isospin breaking effects in the K — > 7r7r(7) system 
of [25] , the differences with the isospin conserving case are 
rather small. In addition we have studied the dependences 
on the subtraction scale /i, where the various assumptions 
are made. Our full amplitudes are /i independent as they 
should. 

The fits show a similar quality to the ones performed 
earlier. The main differences are that the experimental 
Dalitz parameters have changed in K'^ Tr~^iT~^Tr^ and 
are now in better agreement with the ChPT fits. This is 
purely experimental and has nothing to do with the inclu- 
sion of isospin breaking effects. The total is somewhat 
worse because the experimental errors on various partial 
widths have been reduced. 

We also checked how well a few models of the NLO 
weak low energy constants work. The dominance by vec- 
tors and the weak deformation model gave a rather bad 
fit. The factorization model gave a somewhat better fit 
when an extra parameter, an overall scale factor, was al- 
lowed. The quality of this compared to the optimal ChPT 
fits can be judged from the tables giving the best fit values 
for the experimental quantities directly. 
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